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Abstract — Vertical profiles of infrared cirrus extinction have been derived from tropical and 
subtropical upper tropospheric solar occultation spectra. The measurements were recorded by 
the Atmospheric Trace Molecule Spectroscopy (ATMOS) Fourier transform spectrometer 
during the Atmospheric Laboratory for Applications and Sciences (ATLAS) 3 shuttle flight in 
November 1994. The presence of large numbers of small ice crystals is inferred from the 
appearance of broad extinction features in the 8-12 pm region. These features were observed 
near the tropopause and at lower altitudes. Vertical profiles of the ice extinction (km -1 ) in 
microwindows at 831, 957, and 1204 cm" 1 have been retrieved from the spectra and analyzed 
with a model for randomly oriented spheroidal ice crystals. An area-equivalent spherical radius 
of 6 pm is estimated from the smallest ice crystals observed in the 8 -12 pm region. Direct 
penetration of clouds into the lower stratosphere is inferred from observations of cloud 
extinction extending from the upper troposphere to 50 mbar (20 km altitude). Cloud extinction 
between 3 and 5 pm shows very little wavelength dependence, at least for the cases observed by 
the ATMOS instrument in the tropics and subtropics during ATLAS 3. Published by Elsevier 
Science Ltd. 


1. INTRODUCTION 

The frequent occurrence of subvisible cirrus clouds in the tropical upper troposphere 1-4 is of interest 
because of the potential impact of these clouds on the Earth’s radiation budget. Although they only 
slightly increases the albedo, subvisual cirrus can cause a net warming on the order of 1 Wm -2 . 4 
This effect is small if considered only at a single time or location, but it is important because of the 
persistence of these clouds. Thin cirrus IR extinction in the tropopause region may increase the local 
heating rate by a few Kelvin per day, 5,6 resulting in local radiative energy absorption and vertical 
motion that may be in part responsible for lower stratospheric circulation. 6 Furthermore, the 
formation of these clouds may play a role in the freeze drying of air as it is lofted into the tropical 
lower stratosphere. 7 Efficient reactions of C10N0 2 with HC1 and H 2 0, and HOC1 with HC1, are 
likely to take place on the surfaces of cirrus clouds 8 and persistently perturb chlorine chemistry near 
the tropical tropopause, though the changes predicted are small in absolute terms because of the low 
inorganic chlorine abundance. 9 
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In situ observational information on the microphysical properties of ice crystals in tropical upper 
tropospheric cirrus is very limited because these clouds are relatively inaccessible, and it is difficult to 
detect and reliably measure the sizes and shapes of small ice crystals (< 25 pm) with probes. Despite 
considerable scatter, the published in situ measurements" * 14 show a general decrease in size with 
decreasing temperature, with a large number of small crystals observed in the cold tops of tropical 
cumulonimbus systems. Small ice crystals sampled near Kwajalein, Marshall Island, and near Nadi, 
Fiji, made up more than half the mass and caused more th in half the extinction in the upper, colder 
parts of the cirrus. 11 Knollenberg et al 12 noted that the ice crystals in tropical anvils associated with 
tropical cyclones over Micronesia were much smaller than those observed in tropical anvils at 
slightly warmer temperature over Panama. 14 Mean miss equivalent sizes of 20-40 pm were 
measured near the cold cloud tops over Micronesia. 12 Knollenberg et al 12 further noted that the ice 
crystals they observed at very low temperatures were dense and quasi-spherical or of low aspect 
ratio, whereas Hemysfield and McFarquhar 11 noted pr marily columns at temperatures below 
— 50 C. The variety may be due to both natural variabilit y and differences in the sensitivities of the 
instrumentation. 

The ice density in cirrus near the tropical tropopause is over 50 times the water vapour density. 12 
Lidar observations have shown that a narrow (< 1 km ve'tical thickness) subvisible cirrus layer is 
persistent near the tropical tropopause with backtrajectory calculations indicating that this layer 
originates in convective cloud systems. 7 Passive spaceborne instruments have also been used to 
study tropical high altitude cirrus. Microphysical properties have been inferred from observations 
by nadir-viewing infrared imagers 15 * 17 and the Infrared Interferometer Spectrometer (IRIS) 
flown on the Nimbus 4 meteorological satellite in 1970. 18 Nadir-viewing instruments provide 
good spatial and temporal sampling, but very limited vertical information. Limb-viewing instru- 
ments have the advantages of high vertical resolution and sensitivity, but with poorer spatial 
and temporal sampling. To date, the only limb observations analyzed to derive the frequency 
of occurrence and microphysical properties of high altitude thin cirrus are the visible and near 
infrared solar occultation observations recorded by the Sti atospheric Aerosol and Gas Experiment 
(SAGE) II instrument. 2 4 Assuming spherical particles and Mie theory, Wang et al 2 concluded that 
the SAGE II measurements contain information on particles with radii between about 0.1 and 
1.0 pm. The SAGE II tropical high cloud data reveal varying ratios between the extinction at 0.525 
and 1.02 pm, which have been interpreted as indicating ice particle growth during cirrus cloud 
evolution. 

The purpose of this paper is to report an initial study of high spectral and high vertical resolution 
infrared solar occultation spectra of tropical and northerr hemisphere subtropical ice clouds. The 
measurements were obtained from low Earth orbit with tf e Atmospheric Trace Molecule Spectro- 
scopy (ATMOS) Fourier transform spectrometer (FTS) daring its final shuttle flight in November 
1994. We focus here on the analysis of spectra recorded w th filters covering 625-1400 cm * 1 (filter 
12) and 1550-3420 cm* 1 (filter 3). An accompanying paper 19 presents an analysis of the same data 
set for the volume mixing ratio profiles of several relativdy long-lived trace gases. In this work, 
a brief description of the tropical and subtropical FTS measurements is presented in Sec. 2. In 
Sec. 2.1, filter 12 high and low Sun spectra are illustrated foi clear sky conditions and an observation 
through a thin cirrus cloud with a high density of small ice crystals. The method used to retrieve 
vertical cirrus extinction profiles at selected narrow spectial regions is described and examples of 
multi-wavelength vertical cirrus extinction profiles are p-esented. Section 2.2 presents a similar 
discussion and results from the filter 3 observations. In Sec 3, we report the comparison of selected 
cirrus extinction profiles with calculated values generated w ith a model developed for simulating the 
infrared radiative properties of spheroidal and hexagonal ice crystals. 17 The measurements are 
compared with model results calculated for spherical particles and spheroids. Finally, the key 
conclusions and unresolved questions are summarized in sec. 4. 


2. ATMOS/ATLAS 3 TROPICAL AND SUBTt OP1CAL OBSERVATIONS 

The ATMOS FTS recorded broadband 0.01 cm -1 resolution mid-IR solar occultation spectra 
during 4 Shuttle flights. 20,21 Two-sided interferograms with a maximum optical path difference of 
47 cm were obtained at 2.2 s intervals during sunrises or s insets and Fourier transformed to yield 



ATMOS/ATLAS 3 infrared profile measurements of clouds 


905 


Table 1. ATMOS/ATLAS 3 occultations analyzed in this study 


Event 

Date 
1994 UT 

Lai. 
deg N 

Long, 
deg E 

Minimum tangent 
altitude (km) 

Sky 

conditions 

Trop. height 
(km) 

Filter 3 (1550- 3420 cm 







SS70 

November 9 

31.4 

94.0 

9 

CS 

14.9 

SS73 

November 9 

30.5 

26.5 

5 

cs 

13.9 

SS76 

November 10 

26.5 

184.3 

5 

CS 

17.2 

SS79 

November 10 

24.7 

117.0 

7 

LAC 

16.1 

SS85 

November 10 

22.0 

342.2 

11 

LAC 

16.0 

SS88 

November 10 

19.9 

275.0 

13 

LAC 

16.3 

SS91 

November 1 1 

18.0 

207.7 

16 

HAC 

16.8 

SS94 

November 1 1 

15.9 

140.4 

8 

LAC 

16.6 

SS97 

November 1 1 

14.0 

73.1 

18 

HAC 

16.6 

SSI 00 

November 1 i 

11.3 

343.5 

< 7 

LC 

16.6 

SS103 

November 1 1 

8.6 

276.4 

18 

HAC 

16.5 

SS106 

November 12 

6.6 

209.1 

13 

LAC 

16.7 

Filter 12 (625-1400 cm 

l ) 






SS71 

November 9 

31.1 

71.6 

< 5 

CS 

16.0 

SS74 

November 9 

30.2 

4.0 

< 5 

CS 

13.6 

SS80 

November 10 

23.9 

94.6 

15 

HAC 

15.7 

SS83 

November 10 

22.7 

27.2 

3 

CS 

16.2 

SS86 

November 10 

21.4 

319.7 

< 5 

CS 

16.. 2 

SS89 

November 1 1 

19.1 

252.6 

16 

HAC 

16.4 

SS92 

November 1 1 

17.8 

185.1 

17 

HAC 

16.8 

SS95 

November 11 

15.2 

118.0 

6 

LAC 

16.6 

SS98 

November 1 1 

13.0 

50.8 

12 

HAC 

16.4 

SSI 04 

November 12 

8.0 

253.9 

< 5 

LC 

16.6 

SS107 

November 12 

5.3 

186.8 

14 

HAC 

16.8 


Note : Latitude (Lat.) and longitude (Long.) at the tangent point for the spectrum nearest to a tangent altitude of 15 km. 
SS-sunset. Qualitative classification of upper tropospheric measurements: CS-clear sky; HAC-high altitude cirrus (maximum 
altitude > 15 km); LAC-low altitude cirrus (maximum altitude < 15 km); LC-layered cirrus. Tropopause (Trop.) heights are 
from NMC data 4,5 . 


series of exoatmospheric (pure solar) and low Sun (atmospheric) spectra below a tangent altitude of 
^150 km. A two-axis servo-controlled suntracker was used to position the sensor field of view on 
the Sun during the 4 min of observations during each occultation event. Overlapping optical 
bandpass filters limited the total spectral coverage to 625 to 4750 cm - V The results reported in this 
work were recorded during the Atmospheric Laboratory for Applications and Science (ATLAS) 
3 mission (3-12 November 1994) 22 w'ith ATMOS filter 12 covering 625-1400 cm" 1 at a signal to 
root-mean-square (rms) noise ratio of ~ 100. Each low Sun spectrum was divided by an average of 
exoatmospheric spectra from the same occultation to remove solar features, lines of H 2 0 and COi 
from residual air within the interferometer, and the wavelength dependence of the instrument-filter- 
detector system response. 

ATMOS FTS atmospheric observations during the ATLAS 3 mission were recorded during 107 
sunsets with tangent point latitudes between 5°N and 49 N and 94 sunrises with tangent point 
latitudes of 64°S to 72°S. 23 Tropical and subtropical measurements, the focus of this paper, were 
obtained late in the mission at relatively high beta angles. The refracted tangent height spacings of 
successive late mission observations (~1 km near the tropopause) were one-half to one-third the 
typical ATMOS value. The effective vertical resolution of the observations was limited principally 
by the finite scan time and the circular instrumental field of view on the Sun: 1.0 mrad in diameter 
(2.0 km at the tangent point) for filter 3 observations and 2.0 mrad in diameter (4.0 km at the tangent 
point) for filter 12 observations. The range in tangent point latitude and longitude sampled during 
a tropical occultation varied significantly during the event but did not exceed 2 between tangent 
heights of 15 and 30 km. Table 1 provides a list of the tropical and subtropical occultations analyzed 
in this study. Papers by Norton and Rinsland 24 and Abrams et al 25 provide a description of the 
ATMOS data-processing methods which have been applied primarily to retrieve stratospheric and 
mesospheric molecular volume mixing ratio and pressure-temperature profiles from microwindows 
in the spectra. The accompanying paper 19 describes the method used to determine tangent pressures 
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and volume mixing ratio profiles of several molecules from the upper tropospheric spectra analyzed 
in this work. 

2.1. Filter 12 spectra 

Figure 1 shows examples of the 8-12 pm window region in tropical spectra recorded with filter 12. 
Each spectrum has been normalized to the peak signal in the interval. The upper panel presents 

1.2 - — — — 



750 850 950 10 50 1150 1250 

WAVENUMBER (cm" 1 ) 

I 1 



750 850 950 10 50 11 50 1 250 

WAVENUMBER (cm"') 

Fig. I Examples of the 8-12 pm region in ATMOS tropical filte * 12 spectra. Each spectrum is normalized 
to the highest signal in the interval. The upper spectrum in the to] panel is offset vertically by 0.2 for clarity. 
Upper panel, scan 1: high Sun spectrum from SS104 recorded < t a tangent pressure of 2.83 x 10 4 mbar 
(100.0 km tangent height). Upper panel, scan 2: upper troposph< ric spectrum from SSI 04 recorded under 
nearly clear sky conditions. The pressure, attitude, and temperature at the tangent point were 104.0 mb, 
16.3 km, and 198.2 K, respectively. Locations of prominent mole< ular absorption features are marked with 
letters. See text for identifications. Lower panel, scan 3: SSI 07 spectrum recorded through a tropical cirrus 
(ice) cloud. The pressure, altitude, and temperature at the tangent point were 103.3 mb, 16.3 km, and 

197.4 K, respectively. 
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observations recorded during occultation SSI 04 (SS-sunset). The tangent pressures (heights) of the 
scans are 2.83 x 10 4 mbar (100.0 km) and 104.0 mbar (16.3 km), respectively. No significant absorp- 
tion features appear in the upper scan, which illustrates the high signal-to-noise of the observations 
and the success achieved in removing instrumental and solar features from this spectrum. The low 
Sun spectrum was recorded under nearly clear sky conditions. Nine prominent molecular absorp- 
tion features are marked above the spectrum: (A) identifies the 794 cm" 1 v 3 band of CC1 4 ; (B) the 
809 cm 1 v 7 band of CHC1F 2 (HCFC-22); (C) the 846 cm 1 v 4 band of CC1 3 F (CFC-1 1); (D) the 
overlapping 879 cm' 1 v 5 and 897 cm" 1 2\) g bands of HN0 3 ; (E) the 923 cm” 1 v 5 band of CC1 2 F 2 
(CFC-1 2); (F) the 960cm' 1 (000 1 1 M 1 000 1 ) band of C0 2 , (G) the 1042 cm' 1 v 3 band of 0 3 ; 
(H) the 1102 cm' 1 iq band of CC1 2 F 2 ; and (I) the 1161 cm' 1 v 8 band of CC1 2 F 2 . This spectral 
region also contains absorption by weaker bands of the molecules noted above plus significant 
absorptions by additional molecules, particularly N 2 0, CH 4 , H 2 0, and SF 6 . 

The lower panel of Fig. 1 shows a normalized spectrum recorded during occultation SS107. The 
pressure and temperature at the tangent point are nearly the same as the values for the low Sun 
spectrum from SS104. The molecular features in the SS104 and SS107 low Sun spectra have very 
similar absorption depths. However, the molecular bands in SS107 are superimposed on additional 
broad extinction features. A distinct feature with a maximum is observed at 955 cm' 1 . This feature 
has been observed previously in ground-based, 26 nadir-view aircraft, 27 ' 29 and IRIS-D satellite 
nadir-view emission FTS spectra in equatorial 1 and polar regions. 30 The peaked spectral shape has 
been reported previously for this feature only for nadir observations in polar regions during the 
winter and the early spring. 30 Extinction by this feature increases again above 960 cm' 1 , but the 
wavelength dependence is difficult to quantify between 960 and 1080 cm" 1 because of strong 
absorption by the 0 3 v 3 band. The normalized signal between 1080 and 1250 cm' 1 is also 
significantly lower in SS107 than in SSI 04. The additional extinction in this region is overlapped by 
moderately strong molecular absorption by several molecules. 

As noted in previous studies, 26 ' 29 the spectral properties of the additional, broad extinction 
observed in SSI 07 but not in SSI 04 can be attributed to small water ice and/or liquid water particles 
in optically thin cirrus clouds. Ice crystals can be expected at the low temperatures near the tangent 
point of the SS107 observations. This extinction cannot be explained by the presence of sulfate 
aerosol (H 2 S0 4 + H 2 0) droplets, which have a much different wavelength dependence in the 
8-12 pm region. 31 Although HN0 3 is soluble in water and has been predicted to condense with 
water as nitric acid trihydrate (NAT) at temperatures slightly above the frost point, calculations of 
aerosol composition made with an equilibrium model indicate that HN0 3 plays only a minor role in 
the crystal nucleation in the tropical upper troposphere because of the low concentration of HN0 3 
vapour. 5 Furthermore, as demonstrated by a recent reanalysis 32 of the infrared spectra of type I 
polar stratospheric clouds ( PSCs) observed over Antarctica in September 1987, 33 the role of HN0 3 
in low temperature cloud formation is unclear; the observed infrared spectral features do not match 
synthetic spectra computed for NAT and other suspected compositions. A more complex, but 
unknown composition is likely 32 for the aerosols in these low temperature, nitrogen-enriched 
clouds. 34 

Window measurements at 8, 1 1, and 12 pm are useful for distinguishing between ice water and 
liquid water clouds 29 and quantifying the sizes and shapes of small ice crystals in cirrus. 15 " 1 7,26-29 
Therefore, cloud extinction was quantified from measurements in narrow window regions near these 
three wavelengths. The signals from the full set of filter 12 tropical observations were first plotted 
versus tangent pressure. The highest signals at the three wavelengths form a compact relation vs 
tangent pressure which was least-squares fitted and assumed to represent clear sky conditions. The 
signal defined by this curve varies with tangent pressure because of the nonlinear response of the 
ATMOS HgCdTe photoconductive detector at the high incident solar flux. 35 About half of the 
observations below the tropopause have signals lower than the reference curve. We assume that the 
low signals are the result of cloud extinction along the measurement line of sight. 

Figure 2 shows the residuals calculated from the difference between the clear sky reference curves 
and measurements from occultations SS71, SS74, and SS104. The residuals for 830.75-831.00 cm ' 1 
are presented in the left panel. Hereafter, this interval is called the 831 cm' 1 window. Laboratory 
measurements between 100 and 266 K indicate that the imaginary index of refraction of hexagonal 
ice has a peak near this wavelength. 36 ' 40 Line-by-line radiative transfer calculations with a realistic 
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Fig. 2. Signals measured in the 831, 957, and 1204 cm -1 windows divided by the corresponding clear sky 
reference curve. The observations were recorded under nearly optically clear sky conditions during 
occupations SS71, SS74, and SS104. Arrow heads mark two low values from optical thin cirrus measured in 
the tropopause region. Error bars show standard deviations of the values. 


set of gas mixing ratio profiles and the molecular parameters from the 1995 ATMOS line list 41 and 
the 1996 HITRAN compilation 42 (lines and absorption cross sections included from both databases) 
predict a transmittance of 0.98 for limb paths above a tangent pressure of 150 mbar (14 km altitude). 
The measurements marked with arrowheads show the extinction by an optically thin cirrus layer in 
SS104 with a vertical thickness of ~1 km near the tropical tropopause. 5,6 As shown, this same 
extinction feature appears in the residuals for 956.70^-957.00 cm - 1 (957 cm -1 window) and 
1204.40-1204.55 cm -1 (1204 cm -1 window). Radiative transfer calculations with known molecular 
lines, absorption cross sections, and a reference set of vertical gas concentration profiles predict 
a transmittance of 0.99 and 0.97 above 1 50 mbar (14 km abitude) for these two windows, respective- 
ly. Comparison of the clear sky signal profiles from the 204 cm -1 window with those from the 
windows at 831 and 957 cm -1 confirms that extinction ii the 1204 cm -1 window is greater than 
corresponding values in the other two windows. However, the measured signals at 1204 cm -1 are 
lower than those at 831 and 957 cm -1 by ~13% at 150 mbar, which suggests the presence of 
significant, unmodeled upper tropospheric extinction in :he 1204cm -1 region. We assume that 
molecular and other non-cloud sources of extinction in all three windows have been removed to first 
order by dividing the signal measured in the window by th^ corresponding value from the clear sky 
vs pressure reference curve. Based on the residuals in Fig. 2, the 1 -sigma precisions in defining the 
clear sky reference curves at 831, 957, and 1204 cm -1 in the troposphere and lower stratosphere are 
estimated as 3, 3, and 5%, respectively. A major sulfate aerosol band has its peak extinction near 
1190cm -1 . 31 Hence, significant aerosol extinction will oc;ur in the 1204 cm -1 window over limb 
paths with high aerosol loadings. 

Cloud transmittances, calculated by dividing the measured signals by the clear sky reference 
curve, were converted to optical depths and then analyzed to derive vertical profiles of cloud 
extinction (km -1 ) at the wavelengths of the three windows. An onion-peeling inversion procedure 
was adopted assuming a spherically symmetric atmosp tere. Hence, cloud extinction over the 
200 km observation path is assumed to be only a function of the pressure altitude of the layer. 
Horizontal variations across the width of the field of view md nonspherically symmetric extinction 
along the solar occultation ray path cannot be determined from the observations. 
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Fig. 3. T op panel: Examples of cloud transmittance profiles at 83 1 cm ~ 1 calculated form the average signal 
at 831 cm 1 divided by the corresponding value from the clear sky reference curve. Precisions are ±5%, 
1 <j. Middle panel: Profiles of EXT831/EXT957 for the same occultations. Bottom panel: Profiles of 
EXT831 /EXT 1204 for the same occultations. Error bars are estimated l<r uncertainties. 


The top panel of Fig. 3 presents a plot of the 831 cm" 1 cloud transmission profiles from four 
occultations. The middle panel shows the profile of EXT(831)/EXT(957), the ratio of the cloud 
extinction coefficient (km -1 ) at 831 cm -1 to the cloud extinction coefficient (km -1 ) at 957 cm" 1 , 
retrieved from these observations. Values are plotted only for cloud transmissions at 831 cm 1 
between 0.05 and 0.95. Hence, we exclude observations of weak and strong cloud extinction, which 
yield larger measurement uncertainties. The high extinction measurements are also avoided because 
they include a component due to self-emission by the cloud. Ratios of EXT(831)/EXT(957) larger 
than unity are characteristic of extinction by small ice crystals. 15 “ 17 ’ 26 " 30 The SS107 observations 
indicate that small ice crystals were present in large numbers near the tropopause. During SS98, thin 
cirrus extinction was measured between 80 and 200 mbar with an apparent vertical distribution of 
small crystals that peaked in abundance near a tangent pressure of 150 mbar (14 km altitude). 
Extinction ratios near unity were measured for SS89 and SS95, suggesting that small crystals did not 
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dominate the ice particle populations measured near the tropopause and at lower altitudes, 
respectively. The bottom panel of Fig. 3 presents measurements of EXT(831)/EXT(1204h the ratio of 
the cloud extinction coefficient at 831 cm 1 to the corres >onding value at 1204 cm" 1 . The differ- 
ences among the measurements are small and marginally s gnificant, given the estimated error bars. 
The results presented in Fig. 3 illustrate the variations that occur in tropical and subtropical upper 
tropospheric ice cloud infrared optical properties and the advantages of high vertical and high 
spectral resolution for passive remote sensing cloud studies. 



1600 1900 2200 2500 2800 31 00 3400 

WAVENUMBER 
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1600 1900 2200 2500 2800 3100 3400 
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Fig- 4. Examples of the 3 6 jam region in ATMOS tropical filler 3 spectra. Upper panel, scan 1: high Sun 
spectrum from SS100 recorded at a tangent pressure of 1.31 x U)“ 5 nib (126.7 km tangent height), offset 
vertically for clarity. Upper panel, scan 2: SS100 spectrum rccor Jed near the tropopause under clear sky 
conditions. The pressure, altitude, and temperature at the tar gent point were 76.1 mb, 18.1 km, and 
199.9 K, respectively. Prominent molecular absorption feature » are marked with letters. See text for 
identifications. Lower panel, scan 3: SS97 spectrum recorded in i he lower stratosphere through a tropical 
cirrus (ice) cloud. The pressure, altitude, and temperature at the angent point were 68.1 mb, 18.8 km, and 

204.1 K, respectively. 
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2.2. Filter 3 spectra 

Figure 4 presents plots of high and low Sun ATLAS 3 tropical spectra recorded with ATMOS 
filter 3. The measurements in the upper panel were obtained during occultation SS10O. Spectrum 
1 was recorded at tangent pressure (height) of 1.31 x 10“ 5 mb (126.7 km). The only significant feature 
is the intense v 3 band of C0 2 at 2349 cm -1 , labeled D. Spectrum 2 was recorded under clear sky 
conditions at a tangent pressure (altitude) of 76.1 mbar (18.1 km). The identifications and posi- 
tions of the other prominent molecular features are: (A) the wing of the intense v 2 band of 
H 2 0 at 1594 cm -1 ; (B) the (11 102H00001) band of C0 2 at 1932 cm" 1 ; (C) the overlapping 
( 1 1 101 ) --(00001 ) band of C0 2 at 2076cm" 1 , the v x + v 3 band of 0 3 at 2111 cm" 1 , the v 3 band of 
N 2 0 at 2224cm" 1 , and the (1-0) band of CO at 2143 cm" 1 ; (E) the Vj + 2v§ band of N 2 0 at 
2475 cm" l ; (F) the 2v { band of N 2 0 at 2563 cm" 1 , (G) the overlapping + v 2 + v 3 band of 0 3 at 
2785 cm" 1 and the v 2 + v 4 band of CH 4 at 2838 cm" 1 ; (H) the v 3 band of CH 4 at 3019 cm -1 ; and 
(I) the (21 102)— (00001) band of C0 2 at 3339 cm" 1 . 




SIGNAL( 1963 cm _1 )/REFERENCE 

Fig. 5. Upper panel: Residuals from least-squares fits to measurements of average signal in the 1963 cm “ 1 
window divided by the corresponding values from the empirically-determined reference curve. The 
measurements were recorded under optically clear sky conditions during occultations SS70, SS73, and 
SS76. Error bars show standard deviations of the values. Lower panel: Profiles of cloud transmission in the 
1963 cm -1 window from spectra recorded during occultations SS97, SS100, and SS106. The dashed line 
shows the approximate altitude of the tropopause (TROP). Precisions are ±5%, la. 
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Fig. 6. Upper panel: Ratios of cloud transmittances measured in the 3228 and 1963 cm 1 windows plotted 
versus cloud transmittance measured in the 1963 cm” 1 window. Lower panel: Ratios of cloud transmittan- 
ces measured in the 2661 and 1963 cm" 1 windows plotted versus the cloud transmittance measured in the 
1963 cm” 1 window. Tropical measurements obtained between tangent pressures of 50 and 200 mb 
(12 20 km altitude) have been included. Pre visions are ±5%. \a. 


Spectrum 3 in the lower panel of Fig. 4 was recorded through a tropical ice cloud during 
occultation SS97 at a tangent pressure of 68.1 mbar (H.8km tangent height). Despite strong 
extinction (see Fig. 5), which decreased the signal-to-noise c f the spectrum, the normalized spectrum 
is very similar to the normalized low Sun, clear sky spect urn in the upper panel. The absorption 
depths of the molecular bands are nearly the same, and the upper envelopes of the two spectra 
exhibit only minor differences in their wavelength dependence except for several positive “spikes” in 
the SS97 spectrum. These sharp features are solar lines, which appear as false emission features due 
to suntracker field-of-view motion in response to changing ;loud extinction on the lower limb of the 
solar disk. Spectral intervals containing these false emissio i features were avoided in analyzing the 
ATMOS measurements for cloud extinction. 

Cloud transmittance profiles were measured at 1963. 55-1963.40 cm 1 (1963 cm 1 window), 
2660.80 2661.20 cm 1 (2661 cm 1 window), and 3228.20- 1229.00 cm ~ 1 (3228 cm 1 window) with 
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Fig. 7. Mass extinction cross sections calculated with Mie theory and the hexagonal ice crystals optical 
constants of Toon et al. 39 The assumed particle sizes are indicated in the legend. 


the method described in Sec. 2.1. Line-by-line calculations of molecular absorption with a realistic 
set of volume mixing ratio profiles predict clear sky transmittances >0.99 in these windows above 
150 mbar (14 km). The upper panel of Fig. 5 presents residuals of least-squares fits to 1963 cm 1 
window measurements from the clear sky occultations SS70, SS73, and SS76. Based on these results, 
we estimate a precision of ±3%, ler, for the clear sky vs pressure reference curves in the upper 
troposphere and lower stratosphere. It was necessary to select 26 'N-31 'N latitude occultations for 
defining the filter 3 clear sky reference curves because of detectable upper tropospheric and/or lower 
stratospheric cloud attenuation at some altitudes in all of the tropical occultations. The lower panel 
of Fig. 5 presents the 1 963 cm “ 1 window cloud transmittance profiles derived for occultations SS97, 
SS100, and SSI 06. The SS97 measurements show cloud extinction extending across the tropopause 
into the lower stratosphere. Hence, these observations represent an example of direct convective 
transport of upper tropospheric air into the lower stratosphere. The SS100 measurements show 
extinction by layered clouds at lower altitudes. The cloud transmittance decreased sharply to 0.2 at 
200 mbar (12.2 km) followed by a recovery to 0.9 at 230 mb (1 1.1 km). Clear sky conditions occurred 
above 110 mbar (15.9 km) during SS106 with a sharp decrease to a cloud transmission <0.1 below 
160 mbar (13.3 km). 

Figure 6 presents cloud transmittances at 2661 and 3228 cm -1 , normalized to the correspond- 
ing value at 1963 cm -1 , plotted vs cloud transmittance at 1963 cm -1 . Tropical measurements 
between tangent pressures of 50 and 200 mb (12-20 km altitude) have been included. As 
demonstrated by this figure, upper tropospheric cloud extinction shows very little wavelength 
dependence in the 3-5 pm region, at least for the cases observed by ATMOS in the tropics 
and subtropics during the ATLAS 3 mission. In Fig. 7, Mie theory calculations of the extinction 
cross section for hexagonal ice are presented for the 800-4000 cm - 1 region and particles w F ith radii 
of 0.75, 3, 6, and 12 pm. The strong 3 pm extinction feature predicted for a particle radius of 0.75 pm 
was not detected in any of the ATMOS upper tropospheric filter 3 tropical and subtropical 
observations. 
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3. ICF CRYSTAL LXTINCTION CALCULATIONS 

Monochromatic optical depth ratios have been calculated for randomly oriented spheroidal ice 
crystals on the basis of the model of Takano et al. 1 Following the nomenclature used in Ref. 17, the 
equivalent spherical radius of a spheroid may be expressed by 

ct* = it K [v/i/cos 2 C + v 2 sin ’ 0 1/2 ] kl * ( 1 ) 

where C is the angle between the incident direction and the rotation axis of the spheroid; a s is the 
radius of the spheroid, i.e. the length of the short (long) se niaxis for prolate (oblate) spheroids; and 
r is the axial ratio (> 1 for prolate spheroids, < 1 for oblate spheroids). The extinction cross section 
of a spheroid may be expressed by 


CAO = Qsph^x [(cos 2 C + v 2 iin 2 C) 1/2 ] k2 , 


( 2 ) 


where 0 sph is the optical efficiency factor from Mie theory for a sphere of radius a s . The empirical 
constants k 1 and k2 in Eqs. (1) and (2) were determined as described in Ref. 17, where their numerical 
values are given. Finally, the extinction between cross section for randomly oriented spheroids is 
given by 





C e (C) sin £ do 


(3) 


The Takano et al 17 model may be used as a first approximation for spheroids whose area-equivalent 
size parameter, x^ = 2nrj/. is less than -30, with a e and v related to r g , the radius of the 
area-equivalent sphere, as given in Ref. 43. Area-equivaler t spheres are used because they produce 
the same optical depth as randomly oriented nonspherical particles. 17 

Calculations were performed with tw f o sets of water ice optical constants as provided on the 
HITRAN 1996 CD ROM. 42 Imaginary indices of refraction in the compilation of Warren 38 are 
based on spectral reflectance measurements of polycrystalline ice samples formed from the freezing 
of w ater at 266 K 36 with corresponding real indices of refraction calculated from a Kramers-Kronig 
transform based on imaginary indices of refraction measured over a broad frequency range. The 
163 K hexagonal ice optical constants of Toon et al 39 were determined from thin film transmittances 
and an iterative Kramers Kronig technique. Differences among the reported infrared optical 
constants for H 2 0 ice have been summarized by Toon et al 39 except for the recent polycrystalline 
w ater ice optical constant measurements of Clapp et al 40 (not yet included in HITRAN), which cover 
800-4000 cm 1 at nine temperatures between 130 and 210 K. A strong temperature dependence of 
the IR optical constants is indicated by the Clapp et al 40 datasets near the 3 and 12 pm ice bands. 
The following behavior has been noted for 800-1200 era" l , our primary region of interest, by 
comparing measured optical constants for 130 to 266 K, 36 38 40 As the temperature is lowered over 
this range, the imaginary index of refraction increases below- 1000 cm" 1 with a maximum relative 
change of 2.0 at 945 cm" l , near the sharp peak in the ATMOS thin cirrus, small ice crystal spectrum 
in the lower panel of Fig. 1. Above 1000 cm” 1 , where mi limum values of the imaginary index of 
refraction occur due to weak absorption, a more complex temperature dependence has been 
measured with values both higher and lower than at 266 K; at 1050 cm “ 1 , near the minimum of the 
imaginary index of refraction, values at the coldest temperature are more than a factor of two lower 
than at 266 K. Ratios of the real index of refraction reach a maximum of 1.2 at 800 cm” 1 with the 
higher values at 130 K relative to the 266 K measurements possibly due in part to truncation errors 
in the Kramers Kronig integrations. 3 * 08 40 

Figure 8 presents the model-calculated optical depth ra ios plotted as a function of r g for aspect 
ratios v of 0.5 (oblate spheroid), 1.0 (sphere), and 2.0 (prolate spheroid). Differences in the 
831 cm” 7957 cm" 1 optical depth ratio calculated with the two sets of optical constants and the 
values of v in the figure are up to 35% at r g < 1 pm declining to < 14% for 5 < r g < 20 pm. The 
maximum difference in the 831 cm - 7 1204 cm 1 optical depth ratio is 11% with agreement to 5% 
or better for r g > 2 pm. 

Figure 9 presents comparisons of the extinction ratics measured from the filter 12 tropical 
occultations with values calculated with the Warren 38 and ^oon et al 39 optical constants. Error bars 
reflect uncertainties in the cloud transmittances of ±5% at 831 and 957 cm” 1 and ±7% at 
1204 cm 1 due to instrument noise, zero-level offsets in the spectra, and possible deviations of the 
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Fig. 8. Monochromatic optical depth ratios from calculated extinctions at 831, 957 and 1204 cm” 1 , as 
a function of r g , the radius of an area-equivalent sphere for spherical particles (r = 1.0) and spheroids 
(r = 0.5, 2.0). Curves labeled 831/957 refer to ratios of the optical depth at 831 cm -1 to the value at 
957 cm ” l ; similarly, curves labeled 831/1204 refer to ratios of the optical depth at 831 cm” 1 to the value at 
1204cm” 1 . Solid and dashed curves were computed with the hexagonal water ice optical constants of 
Warren 38 and Toon et al, 34 respectively. 


measurements from the clear sky signal vs tangent pressure reference curve assumed for each micro 
window (e.g. Fig. 2). A larger error has been adopted for the cloud transmittance measurements at 
1204 cm" 1 than those at 831 and 957 cm" 1 because of the stronger extinction observed 
at 1204 cm" 1 and the uncertainty in its origin. The observations are in good agreement with the 
calculations, except for the measurements with EXT831/EXT957 > 1.5. The highest EXT831/ 
EXT957 values of 1.8 correspond to a model r g value of pm, with only minor dependence on the 
v parameter. As illustrated by the calculations in Fig. 8, large 831 cm _1 /957cm 1 optical depth 
ratios are related to small crystals. The ATMOS data for the smallest crystals are within the range of 
2-20 pm simulated by Jensen et al 5 for mature clouds near the tropical tropopause. The larger 
measurement-calculation discrepancies for the small crystal observations may reflect the wide range 
of complex shapes often exhibited by ice crystals in cirrus clouds and the limitations inherent in 
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Fig. 9. Comparison between ATMOS ice extinction ratios at 8 31, 957 and 1204 cm * 1 and model-gener- 
ated values obtained with the hexagonal ice optical constants o! Toon et al 3g (upper panel) and Warren 38 
(lower panel). Solid symbols and error bars show the individual measurements and their \a uncertainties. 
Included observations were restricted to tangent pressures (heights) of 70- 150mbar (14-18 km) and ice 
transmittanccs between 0.05 and 0.95 to mimmize the uncertainties. 


assuming a single crystal size rather than a distribution of fizes in the model calculations. Based on 
the Mie theory calculations and the Toon et al 16 ice opt cal constants, we infer from normalized 
extinction kernel calculations 2 that the ATMOS 8-12 pn spectra contain useful information on 
particle radii between about 1 and 30 pm with good sensitivity in the range 2-10 pm, whereas the 
ATMOS 3-5 pm spectra contain useful information on crystals with radii between about 0.3-1. 0 pm 
and good sensitivity for particle radii between 0.4 and 4 y m. 

Mie theory calculations assuming 6 pm for the charactf ristic radius of small ice crystals predict 
a strong spectral feature with extinction that increases by a Factor of 2.8 from 3500 to 3660 cm " 1 (see 
Fig. 7). Unfortunately, atmospheric transmission is low in that region due to strong absorption by 
C0 2 lines and the strong v x band of HNOj, not included in the HITRAN 1996 database. 44 A search 
for this characteristic ice feature in the ATMOS filter 4 spt ctra (3120-4750 cm' 1 ) was inconclusive 
because of the strong molecular interferences. A definitive search will require parameters for the 
HN0 3 v j band. 
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4. SUMMARY AND CONCLUSIONS 

In this paper, we have reported an initial analysis of cloud extinction measurements from 
0.01 cm” 1 resolution, broadband IR solar occultation spectra recorded in the tropical and sub- 
tropical upper troposphere and lower stratosphere by the ATMOS FTS during the ATLAS 3 shuttle 
mission of November 1994. Several occultations show wavelength-dependent ice extinction in the 
8-12 pm region due to large numbers of small ice crystals. None of the occultations observed with 
a filter covering 3-5 pm show the 3 pm ice band, which is a characteristic extinction feature of ice 
particles with radii < 1 pm. Methods have been developed to retrieve cirrus extinction coefficient 
vertical profiles from measurements in narrow spectral intervals with minimal molecular absorption. 
The strength of this approach is that it allows us to take advantage of the high spectral resolution to 
avoid interfering lines and the multiple observations to define clear sky reference curves to 
accurately account for unmodeled extinction. Comparisons of our ATMOS 8-12 pm cirrus extinc- 
tion measurements with calculations for randomly-oriented spheroidal ice crystals indicate an 
area-equivalent spherical radius of about 6 pm for the smallest observed ice crystals. Differences 
between the small particle extinction measurements and calculations performed with the hexagonal 
water ice optical constants of Warren 38 and those of Toon et al 39 imply the need for further study. 
Recent laboratory measurements 40 show a strong temperature dependence for ice optical constants 
in the spectral regions used in this work and many previous remote-sensing cloud 
studies. 1 15-1 7,26-30 Hence, future 8-12 pm cloud investigations need to incorporate temperature- 
dependent optical constants in models of cirrus radiative properties and algorithms used to retrieve 
particle size distributions from atmospheric measurements. Further studies of the ATMOS cloud 
extinction measurements are in progress. 
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